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Abstract 



The normal state of the Holstein model is studied at half-filling in infinite 
dimensions and in the adiabatic regime. The dynamical mean-field equations 
are solved using perturbation expansions around the extremal paths of the 
effective action for the atoms. We find that the Migdal-Eliashberg expansion 
breaks down in the metallic state if the electron-phonon coupling A exceeds 
a value of about 1.3 in spite of the fact that the formal expansion parameter 
Xojq/Ef (ujq is the phonon frequency, Ep the Fermi energy) is much smaller 
than 1. The breakdown is due to the appearance of more than one extremal 
path of the action. We present numerical results which illustrate in detail 
the evolution of the local Green's function, the self-energy and the effective 
atomic potential as a function of A. 
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I. INTRODUCTION 



Studies on the microscopic mechanism of high-T c superconductivity and of the colossal 
magnetoresistance have shown the need for a more detailed investigation of simple models 
for the electron-phonon coupling. One generic model is the Holstein mo delB which con- 
sists of conduction electrons interacting locally with optical phonons. The case of just one 
electron has received a lot of attention in the past. In particular, the formation of small 
and large polarons and their motions have been studied considering both the adiabatic and 
the anti-adiabatic casesEHl. However, the application of the obtained results to high-T c 
superconductors! or to the manganitesll is problematic. In these systems the Fermi energy 
Ep is substantially larger than typical phonon energies ujq. One thus deals in these sys- 
tems with a large Fermi surface where the Migdal-Eliashbergi self-consistent perturbation 
theory, governed by the small parameter Xuq/Ef (A is the electron-phonon coupling con- 
stant), should be more appropriate than the use of canonical transformations typical for the 
one-electron problem. 

From the studies on the single polaron problem it is known that polaron formation, 
i.e., the sudden and large mass enhancement of the electron, takes place around A ~ 1. 
Associated with this is a breakdown of simple perturbation theory due to the appearance 
of a bound state in the electron-phonon system. Whether something similar happens in the 
case of many electrons and a large Fermi surface is still a matter of controversy. On the one 
hand, it has been argued that, due to the Migdal theorem, self-consistent theory works even 
if A is not small. It is claimed that the only condition for the validity of the self-consistent 
perturbation theory is that Xu /E F < 1 holds!. As a result many treatments in the past, 
for instance, for high-T c superconductors, used values for A which are much larger than 10. 
On the other hand, it has been claimed that Migdal-Eliashberg perturbation theory breaks 
down if A exceeds a critical value somewhere between 1 and 2, irrespective of the electronic 
density and the size of the Fermi surfaceB The arguments for this, however, were based 
on the Lang-Firsov transformation^ which is more appropriate in the anti-adiabatic than 
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in the adiabatic limit and they seem to be supported by numerical calculations^. One aim 
of our contribution is to study the possible breakdown of Migdal-Eliashberg perturbation 
theory and polaron formation in the adiabatic regime where the hopping and the Fermi 
energy are much larger than the phonon energies. Another aim is the calculation of phonon 
renormalization effects caused by an arbitrarily strong electron-phonon coupling. 



A suitable method to study the above problem is the dynamical mean-field theoryMEj. 
This method becomes exact in the limit of infinite spatial dimensions. It allows to map the 
lattice problem onto Anderson's single impurity problem embedded in an effective bath cre- 
ated by all the other sites and interacting with the electrons on the impurity in a dynamical, 
time-dependent way. In spite of the substantial simplifications introduced by infinite dimen- 
sions the resulting equations are in general still too complicated to admit exact analytic solu- 
tions. For the Holstein model some results from Monte Carlo simulations are availableEHI. 
Various approximate treatments have also been performed such as iterated perturbation 
theoryS, self-consistent perturbation theory^ and the semiclassical approximation^!!. 

In the following we will formulate the dynamical mean-field equations in terms of path 
integrals. The phonon coordinate will not be treated as a classical variable without dynamics 
as in the semiclassical approach. Instead, we expand the path integrals around the paths 
which extremize the action. For small A there is only one extremal and trivial path and 
Migdal-Eliashberg expansion is valid. Beyond a critical value for A more than one extremal 
paths exist signalizing the breakdown of Migdal-Eliashberg theory though the system is still 
in the metallic state. Our treatment can be applied in the adiabatic regime to all coupling 
strengths A and contains correctly the weak- and the strong-coupling cases. The paper 
is organized as follows: In section II the model and the dynamical mean-field equations 
are formulated. In section III the electrons are approximately integrated out in the path 
integral using a gradient expansion in time. The resulting effective action for the atoms 
is discussed. Approximate expressions for the electron Green's function are obtained in 
section IV. Section V contains numerical results for a Bethe lattice both for the atomic and 
the electronic properties. This section illustrates in detail the evolution of basic quantities 
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of the dynamical mean field theory such as the local Green's function, the electronic self- 
energy, the effective "unperturbed" Green's function and the effective atomic potential as a 
function of A. Finally, section VI contains the conclusions. 



II. MODELS AND DYNAMICAL MEAN-FIELD EQUATIONS 

The Hamiltonian for the Holstein model can be written asQ 

H = ~ £ k 3 {ctc ](J + h.c.) + i E( Wi) 2 + ^ 2 ) + 9 E (!) 

ija i ia 

It describes a system of electrons interacting locally with dispersionless phonons. i,j 
denote the sites of a lattice, u is the frequency of the phonons, tij is the hopping matrix 
element of the electrons, and g the coupling constant between electrons and phonons. 0j is 
the displacement operator of the atom at site i which can be expressed in terms of bosonic 
creation and annihilation operators 

' :(a, + a+). (2) 



c i(T , Ci >a are creation and annihilation operators for electrons at site i and spin direction a 
and n i<T = c} irj Ci a is the density operator of electrons. The partition function of the model 
can be written as a path integral over anticommutating Grassmann numbers c*, c (replacing 
the electronic operators), and phonon fields, 

Z = J D(f)Dc*Dc e s[c '' c ^\ (3) 

where S is the action associated with the Hamiltonian (1). 

In the limit of infinite spatial dimensions, d — > oo, the dynamical mean- field approxi- 
mation for H becomes exact. As a result the many-body problem Eq.(l) can be mapped 
onto an impurity model embedded in a dynamical bath0. Denoting the one-electron Green's 
function associated with the Hamiltonian of Eq.(l) by G(k, iu n ), the self-energy E is defined 
via Dyson's equation 



4 



G-\k,iu n ) = G' »"(k )H ) - S(iw n ). (4) 

, iu; n ) is the free electron Green's function associated with the hopping term in Eq.(l). 
Using the Matsubara formalism uo n are fermionic Matsubara frequencies u n = (2n + l)irT 
where n is an integer and T the temperature. A very important simplification in infinite 
dimensions results from the fact that the self-energy E is independent of the momentum k 
and depends only on the frequency^. The self-energy £ also represents the self-energy of a 
site- independent impurity problem with the effective action S e ff given by 

rP rP 



Seff = I drj dT'J2cUT)G \T-T')c a (r') (5) 

J Q- 



Go describes a dynamical effective field acting on the impurity electrons due to the 
presence of all the other electrons in the crystal. It plays the role of a bare Green's function 
in S e ff and contains all the information on the influence of the surroundings on the impurity 
electrons. (3 in Eq.(5) stands for the inverse temperature 1/T. Denoting the exact one- 
electron Green's function associated with S e jf by Gi oc Dyson's equation of the impurity 
problem reads 

Gj£(iu) n ) = G^iiUn) - E(io; n ). (6) 

Finally, the so far unknown function Go is determined by the condition that G[ oc coincides 
with the local Green's function calculated from the lattice Green's function G(k,iu n ), i.e., 

G loc ^n) = [de , ^ (6) (7) 

J iuj n - e + fjL - S(z^ n ) 

fio is the bare chemical potential and N(e) is the bare density of electron states associated 
with the original hopping term in Eq.(l). In the following we will mainly consider the 
semicircular density of states of the Bethe lattice!! given by 



This lattice has an infinite coordination number z — > oo and a rescaled nearest-neighbor 
hopping term = t/y/z . In this particular case the relation between Gq and Gi oc becomes 

Go 1 (iu n ) = iu n + fiQ - t 2 Gi oc (iu n ). (9) 

Another popular density of states function is that of the Lorentz model 

"< e > = WTW-y (10) 

It describes a model with infinite range hopping terms and a function Go which simply is 
given by 

G (iw n ) = iu n + /i + itsgnu n . (11) 

Go decouples in this case from all the other quantities. This model is therefore especially 
simple allowing often analytic expressions which capture essential features of models with 
more realistic density of states functions. On the other hand it misses the feedback from the 
other quantities and thus is unable to describe, for instance, the Mott-Hubbard transition 
from a metal to an insulator. 

Eqs.(5)-(7) represent a closed set of equations. They have been derived under the as- 
sumption that no long-range order is present which will be assumed in the following. The 
above equations can be solved for large regions in parameter space numerically using Monte 
Carlo or diagonalization methods^. In the following we will use a more analytic method 
which gives directly some insight into the underlying physics. 

III. EFFECTIVE ACTION FOR THE ATOMS 

According to Eq.(5) S e ff is bilinear in the electronic variables. This means that the 
integration over electrons in the partition function Eq.(3) can be carried out yielding the 
following action S p h for the phonon field alone: 

Sph = ~\ [ dT<j>if){~ + uD^r) + Tr log(G 1 - g<j>). (12) 
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Here Tr means the trace over electronic variables, i.e., an integration over r from to (5 
and a sum over the two spin components. The argument of the logarithm in Eq.(12) is 
a nondiagonal matrix consisting of the non-diagonal part Gq 1 and the diagonal part 0. 
Similarly, in u; n -space the first contribution Gq 1 would be diagonal but the second one, 
would become non-diagonal. 

Defining the expectation value of in the usual way by 



< >= J D(j)(f)e^ h /Z, (13) 

with 

Z = J D(j)e s P h , (14) 

there is no reason to expect that this static expectation value is zero even in the half-filled 
case. Indeed, writing down the equation of motion of Heisenberg operators using the original 
Hamiltonian Eq.(l) and considering the static limit we have 

< >= -\ <n> . (15) 

< n > is the average number of electrons per site. It is therefore convenient to split into 
a static and a fluctuation part according to 

=< > +0, (16) 

where < > is defined by the condition < >= 0. This shift can be interpreted as a Hartree 
term T, H = g < > in the self-energy of G . Since this term is frequency independent it 
just renormalizes the bare chemical potential as 

// = fx - S H . (17) 

From now on we will assume that this shift has been performed and that the chemical 
potential in Gq is \i and no longer the bare one /io- 

The main problem left is the evaluation of the Tr log term in Eq.(12). For this we note 
that two independent parameters can be naturally defined in the Holstein model: one is the 



electron-phonon coupling constant A = g 2 /cu^t, and the other one is the adiabatic parameter 
7 = ouo/t. The adiabatic regime is defined by 7 <C 1, corresponding to the situation where 
the electrons move much faster than the phonons. For the calculation of the Tr log term 
we assume that 7 < 1 holds, (f) varies then slowly as function of r compared with the 
variations in the electronic coordinates. An expansion in time gradients of <fi represents thus 
an expansion in the Migdal parameter 7 and we can confine ourselves to the lowest-order 
contributions. 

Expanding the Tr log term in powers of g one obtains 

00 -1 

Tr logiG, 1 - g4>) = Tr logG, 1 - £ -Tr(G g4>) n . (18) 

The first term on the right-hand side of Eq.(18) is independent of and thus will be dropped 
in the following. Writing the conservation of frequencies explicitly by means of a delta- 
function the general term of n-th order in g of the Tr log term becomes 

X n = -2 y —T ]T / dre i ^ + - + ^^(u 1 )G (iy + u 1 )- (19) 

v 

4>(uJ 2 )Gq(v + LOi + u)2)—<i>{u) n )Go(v + Ui + ... + U) n ). 

The prefactor 2 in Eq.(19) is due to the sum over the two spin directions. Due to the delta- 
function there are no restrictions in the sums over uj\...uj n in Eq.(19). For 7 1 the motion 
of the atoms is slow compared to that of the electrons. In leading order in 7 the electrons 
interact thus elastically with the phonons, i.e., they see the momentary value of the phonon 
field as a static field. This means that the phonon frequencies uj\...uj n can be neglected in the 
electron propagators Go. The sums over LO\...uj n can then trivially be carried out in Eq.(19). 
Summing also over n and adding the shifted harmonic potential the total potential V(0(r)) 
in S p h becomes in this local approximation 

Vtitr)) = 4(Hr) + <0» 2 - 2T]T e^ + log(l - gfo)G (iv))- (20) 

In the last term in Eq.(20) we have added an exponential in order to regularize the sum 
over Matsubara frequencies correctly. This regularization affects only the first-order term in 
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g. We found that it is not advisable in general to carry out the frequency sum in Eq.(20) 
directly on the computer but to rewrite it as a contour integral in the complex plane and to 
evaluate it as an integral along the cut on the real frequency axis. One obtains in this way 



2 

2 r+°° 



IT J-oo 



duf(u>) {arctan (^|gg|y) - 7r9(g^G (co + irj) - l)sign$)} . (21) 



Here / denotes the Fermi function, dt and S real and imaginary parts, respectively, and 9 
and sign the theta- and sign-functions, respectively. In the case of the Lorentz model the 
integral in Eq.(21) can be performed analytically and one obtains at zero temperature: 

V{i) = \^o 2 S+{(j>)f - 1/tt[-2# arctan (j^j 

+ t\og (—J——) + n(g4>)sign($) - ngft. (22) 

A conventional way to calculate corrections to the above local terms is to use a gradient 
expansion^. One expands then the propagators in powers of frequencies ui....uj n and rewrites 
these powers in terms of time derivatives acting on the field (p. Such a procedure is, however, 
not possible in our case because of the appearance of nonanalytic terms. To see this let us 
first look at the lowest-order term n = 2: 

2 2 

X 2 = - 9 -Tr{G^G^) =~E 0(n)II(^„)0(-n), (23) 

with 

U(iu n ) = 2 G Q (iu m )G (iuj m + iu n ). (24) 

m 

Using the spectral representation 

G (iu n ) = I de-^- (25) 

and carrying out the integrations one obtains at small frequencies and low temperatures and 
also disregarding a constant term 
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U(iu n ) = 27rp 2 (0)/3| Wri |, (26) 

where p(0) is the spectral function at the chemical potential /i — 0. Considering X 2 as an 
self-energy contribution to the unperturbed phonon propagator in Eq.(12) and performing 
the analytic continuation \uj n \ — > \J —{yj + irf) 2 — > — one recognizes that X 2 describes 
the damping of phonons which is smaller than the phonon frequency u by the Migdal 
paramter 7. Eq.(26) shows, however, that such a damping term produces a nonanalytic 
dependency in S p h on Matsubara frequencies. This rules out a straightforward application 
of the usual gradient expansion. Such an expansion would correspond to a Taylor expansion 
of II in powers of iuj n . Since II has to be even in iu n its linear term would vanish and the 
leading term would be of second order in iu) n in disagreement with the exact result Eq.(26). 
To overcome this problem we will use in the following a modified version of a gradient 
expansion which is able to give the leading non-local corrections to the local term also in 
our case. 

In leading order in 7 or, equivalently, in the local approximation, we assumed that the 
electrons in the general n-th order term X n given in Eq.(19) interact elastically with the 
phonons. The first correction to that obviously is obtained by considering in this diagram 
2 vertices where the electrons are inelastically and n-2 vertices where the electrons are 
elastically scattered. The case of one vertex with inelastic scattering is not possible due to 
the conservation of Matsubara frequencies. The analytic expression for X n is, taking the 
above processes into account, 

n n rP ~ 

*f = --E E T ^e^ + - + ^0(c 1 )Go(^)0(c 2 )G o (^)... 

7~J v 

4>(ui)G (iu + Ui)^(u i+1 )G (iv + u i )...(p(u j )G (iu)^(uj j+1 )G ( y iu)...(j)(uj n )G (iu). (27) 

i and j are the two vertices where inelastic scattering occurs. A factor 1/2 has also been 
introduced in order not to overcount the contributions. However, this factor is compensated 
by the spin factor 2. The superscript 2 at X indicates that all processes are taken into 
account where inelastic scattering takes place at two vertices. 
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The sums over LO\...iOi and Uj + i...u n can be carried out in Eq.(27) yielding 

= -g 2 /nY, T /^rfre i( ^ + ^ )T 0(^)0(cu i )G' (n ~ J ' +l - 1) (^)G' ( ^ l - 1) . (28) 

G is denned by 

G-V) = G o V)-s0(r). (29) 

G^ denotes G, calculated in i-th order in the coupling constant g, i.e., 

G®(iv) = (gfotfG* 1 ^). (30) 

Keeping only the term % — 1 in Eq.(28) is equivalent to omitting the prefactor 1/n. Moreover, 
the sum over j ^ i can easily be carried out yielding 

Xi 2) = -9 2 E /^re^ 1+ ^^0(cu 1 )0(cu 2 )T^G(2z/)G'(2z/ + 2a; 1 ), (31) 

till f tin -,, 



where all terms of order n — 2 have to be taken into account in the product of the two 
Green's functions G. Performing the sum over v in Eq.(31), using the spectral function p of 
G, and also summing over uj 2 we obtain 

*P = -9 2 T [" drdr'(j:^ l{T - T ' ) M)#t)£(t'W(t)], (32) 

with 

F[0]) = -vrp 2 (0). (33) 

p is the spectral function of G which also depends on 0(r). Splitting the sum over oj\ in 
Eq.(32) into positive and negative parts, representing uj 1 as a time derivative which then is 
subjected to a partial integration, and finally performing the sums over u)\ and n we obtain 
for the leading term of the nonlocal part Vnl(<P) of S p h the following expression 

V Kl (i) = f W fWl-^ A,), (34) 
2 Jo 1 — cos{2,ix± (t — t')) or 

with 
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G[4>(r), 0(r')] = F[}(t)] + F[$(t')] - ^— • }(/). (35) 

d(f)(r') 

Vnl vanishes for a constant field so it represents a nonlocal correction to the local potential 
V. Vnl, on the other hand, still depends on two arguments r and r' which would not occur if 
a usual gradient expansion would be applicable. This feature is caused by the nonanalyticity 
of I! in Eq.(26): |a> n | cannot be represented just by one derivative but only by infinite many 
of them causing the appearance of two arguments r and r 1 in Eq.(34). We also would like 
to point out that the expression Eq.(34) is very similar to that derived by Hamann@ for the 
Hubbard model using asymptotic approximations for Go and avoiding a gradient expansion. 



IV. THE LOCAL ELECTRONIC GREEN'S FUNCTION G LO c 

The exact local Green's function Gi oc of the impurity problem can be obtained from the 
relation 

n ( >\ 8lo 9 Z (oc\ 

GUr - T) = 5G Q \r-riy (36) 

Carrying out the functional derivative we obtain 

G loc (r -t')= J DfcGz\n - r 2 ) - gfa x )8(r x - r 2 ))^ T ,e s ^ /Z. (37) 

S p h is given by Eq.(12) where the Tr log term has been approximated by the local and 
nonlocal potentials V and Vjvx. The simplest approximation for the functional integral in 
Eq. (37) is obtained by considering as a classical field and assuming the mass of the atoms 
to be infinite. The kinetic energy of the atoms as well as the nonlocal part in S^h drops then 
out and the functional integral becomes a usual integral. As a result we have 

roo „ „ - 

Giodr -r')= ~ g^e"^ / Z. (38) 

J oo 

Eq.(38) says that G\ oc is to be calculated for a general displacement (p and then the result 
is averaged over all displacements using a Boltzmann factor. In this approximation the 
electrons no longer form a Fermi liquid and there is no quasiparticle peak in the electron 
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Green's function. For a finite mass of the atoms such an approach is valid if T ^> ujq. The 
dynamical mean field approximation on a Bethe lattice amounts then to solve Eqs. (9), (20) 
and (38) self-consistently011. 

A general approximation scheme for path integrals valid especially at low temperatures 
is based on expansions around stationary paths of the action S e ff. These paths are solutions 
of the equation 

= 0. (39) 

S<f>{T) 

We will see that Eq.(39) has in general 3 solutions, two of them being constant in time and 
one depending on time. The surroundings of these solutions will yield the major contribu- 
tions to Gl oc . 

A. Contributions to Gi oc from small fluctuations 

Time-independent solutions of Eq.(39) satisfy also the equation 

*W - < 4o > 

0(f) 

Taking stability considerations into account the relevant solutions of Eq.(40) correspond to 
the minima of V which will be denoted by (j) a , where a counts different minima. We expand 
now the field (f) around each minimum (f) a : 

i = k + <t>' a - ( 41 ) 

Inserting Eq.(41) into S p h we write 

SAk + 4>' a ) = S$ >a + S> phja , (42) 

where S^ a contains the kinetic energy of the atoms and the quadratic terms in <j)' a of V. 
S' ph a is due to the anharmonic part of V around the minimum a and the nonlocal part Vjvl- 
Expanding the inverse in Eq.(37) in powers of <f>' a as well as the exponential in powers of 
S' h a one obtains path integrals over the displacement field containing an arbitrary number of 
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powers in 4>' a . Since the action in these path integrals is S$ a Wick's theorem can be applied. 
Contributions to Gi oc from the minimum a, called G a in the following, are thus obtained in 
terms of usual diagrams due to the electron-phonon and anharmonic interactions. As usual 
it is sufficient to consider self-energy diagrams S Q to G a . G a and S a are related by 

= G ,a ~ (43) 
where Go, a is the shifted unperturbed Green's function defined by 

= Go 1 - 9~K- (44) 

Fig. 1 shows the lowest-order skeleton diagrams for S Q which are of second (diagram 
(a)) and of fourth order (diagram (b)) in g. All uncrossed diagrams are taken into account 
by diagram (a), whereas diagram (b) represents the lowest order vertex correction to (a). 
In conventional treatments of the electron-phonon coupling it is argued that diagram (b) is 
smaller by the factor \uj /t compared to the diagram (a). Whether this is true also in our 
case is not apriori clear because Migdal's theorem relies heavily on properties of momentum 
integrals and phase space arguments for which there is no analogue in our case. 

In order to judge the convergence of the perturbation expansion around (fi a we have 
calculated diagrams (a) and (b). One obtains 

9^ 9^ f 1 
^(iuin) = -- — b(w )G a (iw n + uj ) + 7^ — b(-uj )G a (iuj n - u ) + g 2 \ dep a (e) f (e) : — ^ 

(45) 

where p a {e) is the spectral function of G a and b is the Bose distribution function, 
is given by 

E<M = V (_i)(-+*+2)/2 l_ f dede > de » p a (e)p a (e')p (e") 

4c^7 (-e + iuj n + suj )(-e" + iuj n + tuJo) 

f A(e',s,t) + B(e,e',t) + C(e',e",s) + D(e,e',e") v 

^iuj n + sujq + tuj — e' e — e' + tuj e" + slu — e' e — iuo n + e" — e' ' ' 

with 
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A(e',s,t) = (b(tu ) + f(e'))(b(su; ) + f(e'-tu; )), 



(47) 



B{e,e',t) = (b(tu; ) + /(e'))(/(e) - f(-tcu + e')), 



(48) 



G(e', e", s) = (f(e") - f(e'))(b(su; ) - b{e> - e")), 



(49) 



D(e,e',e") = (/(e") - f(e'))(f(e)+b(e'-e")). 



(50) 



Figs. 2 and 3 present numerical results for the imaginary parts of diagrams (a) (solid lines) 
and (b) (dashed lines) as a function of frequency. In Fig. 2 (3) we used the parameters 
A = 1.0 (1.375), lu = 0.1, and T = 0.01, where all energies are expressed in units of t. 
G a has been calculated using the self-consistently determined function Go described in the 
next section. The potential V develops a double well around A = 1.2 so that Figs. 2 and 
3 correspond to the cases with one and two minima , respectively. If Migdal's theorem 
would hold diagram (b) should be smaller by one order of magnitude compared to diagram 
(a). The Figures show that this is approximately the case which validates the employed 
perturbation expansion. 

The above result is nontrivial, especially, at large couplings when the double well has 
been formed. As illustrated in the next section the main weight of G resides in this case in 
a narrow, central component at zero frequency. The characteristic energy of Go thus may 
be low and comparable with ujq making a simple perturbation theory of S Q in terms of Go, a 
impossible. On the other hand, G a consists at large couplings of a broad peak of width 
~ t around (f> a . The effective expansion parameter for an expansion of X in terms of G a 
can easily be determined for large couplings where one can assume that ImG a {ui) is only 
non-zero for u < —uq. Eqs.(45) and (46) reduce then at T = to 




2u; 



G a {u + ir) + ouq), 



(51) 




G 2 a (uo + ir] + uj )G a (uj + irj + 2uj ). 



(52) 
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The ratio /E^ is thus ~ g 2 joj^t 2 = Xcuo/t which is in agreement with the above numerical 
calculations. Nonvanishing potential minima at 4> a signalize the appearance of nonadiabatic 
effects: at large couplings electrons occupy predominantly ionic configurations with zero or 2 
electrons at the atom. In each case the lattices relaxes to a non- vanishing value similar as in 
the anti-adiabatic case. As shown above, the perturbation expansions around the potential 
minima are, however, still controlled by the Migdal parameter \uo/t. 

According to our previous discussion the nonlocal potential V/vl is smaller by the Migdal 
parameter 7 compared to to the local potential. To make this more quantitative we have 
calculated the imaginary-part of the self-energy of Vnl for a harmonic phonon with frequency 
co>o for a fixed, but general field value <fi. The result is shown in Fig. 4 for the two coupling 
strengths used in Figs. 2 and 3 using again the self-consistent function Go- The position 
of the potential minima are indicated by arrows. One concludes from Fig. 4 that the 
damping of the oscillations around the potential minima is small enough in each case so 
that one may neglect it. Finally, we have neglected anharmonic effects. This is an excellent 
approximation at small and large couplings A but not for intermediate couplings where the 
double well forms. A simple prescription to take anharmonic effects into account which is 
based on the exact atomic eigenf unctions in the potential V will be described in the next 
section. 



B. Instanton contribution to G\ oc 

At half-filling the potential V is symmetric with respect to <fi = 0. If V has two minima 
they are equivalent and the potential barrier between them is finite. The appearance of 
tunnelling processes between the minima manifests itself in a third, time-dependent solution 
of Eq.(39)i. This instanton solution passes from one minimum to the other as function of r 
and is relevant at low temperatures where this passage occurs within the interval [0, f3\. One 
instanton solution creates, however, a whole family of approximate solutions of Eq.(39): 
The passage from one to the other minimum can take place everywhere in the interval 
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[0, /3\; furtheremore, multi-instanton configuration are also approximate solutions of Eq (39) 
and must be taken into account. For each instanton configuration we expand the inverse 
in Eq.(39) in powers of <ft and sum then over all configurations using the rules derived in 
RefS. The first-order term in <p vanishes by symmetry. The second-order self-energy can 
be represented by the diagram a) in Fig. 1 where ther phonon propagator is to be replaced 
by the instanton propagator 

~ 1 - e~ AE/3 2AE 

D ^n) = ^Y^m ■ j$MKEf (53) 

cj)\ is the squared field value of one of the two minima. If V^l can be neglected AE denotes 
the tunnelling frequency of one atom through the potential V between the minima. In 
general, AE contains an additional exponential factor with the action due to Vnl along 
one instanton in the exponent. Dj has the form of a phonon propagator with the phonon 
frequency replaced by AE and an additional prefactor ~ AE/T at temperatures large 
compared to the tunnelling frequency. 

The second-order instanton contribution to the self-energy becomes 

Ej(tuj n ) = 9 2 <P L, _ AE(3 (-b(AE)G (tuj n + AE) + b(-AE)G (iuj n - AE) 

+ - / dep(e)f(e)- . ), (54) 

u J (e - tuj n ) 2 - (AE) 2 s 

where p(e) is the spectral function of Go(e). In the limit AE — > £/ approaches 4>\Go{iuJo). 
Inserting this into Eq.(37) one finds the contribution 

G\ QC = \( : 1 - + : 1 - ), (55) 

to G[ oc . G' loc is due to the integration over the paths along the two minima of V which are 
also included in £/. This contribution has already been taken into account when integrating 
over paths near the first two stationary path at <p a and thus should be subtracted again. 
The pure instanton contribution to G\ oc is thus given in lowest-order perturbation theory by 

Gi = r J - - G' loc . (56) 
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After the analytic continuation the imaginary part of £/ deviates from <ftlGo(u + irf) only 
at small frequencies uj ~ AE where it rapidly decreases and approaches zero at zero tem- 
perature assuming that the characteristic energy scale t* of G is larger than AE. Thus 
the instanton contribution to Gi oc provides the finite value of Gi oc of Fermi liquid theory at 
/i — in form of a Kondo resonance peak with an exponentially small width determined by 



The expansion of the denominator on the right-hand side of Eq.(37) in powers of g<f> 
seems to be problematic for instanton field configurations because in this case g<p can in 
general not be considered as small compared to Gq 1 . This is also true if the perturbation 
series is resummed by means of the self-energy £/ and if the latter is approximated by its 
lowest-order expression in g. In the remainder of this subsection we will show that Gi can 
be obtained in closed form without using an expansion of the denominator in Eq.(37). The 
only assumption is that the instanton gas is dilute, i.e., that self-energy effects for the tunnel 
modes can be neglected. 

Let us first disregard the nonlocal part Vnl in the atomic potential. The Hamiltonian for 
the impurity atom is then hermitean and can be diagonalized yielding two low-lying states 
which are split by tunnelling processes. Disregarding small oscillations around the extremal 
paths the evaluation of the path integral by instanton solutions corresponds to the solution 
of the Schroedinger equation for the two lowest-lying states of the Hamiltonian. Thus the 
instanton part of the electron Green's function, Gj, can also be obtained from an effective 
two-level model with the Hamiltonian 



X 11 and X 00 are Hubbard operators defined as the projection operators |1 >< 1| and 
|0 >< 0|, respectively, where |0 >, |1 > are the two lowest states of the atomic Hamiltonian 
and AE their splitting. The field variable becomes 



AE. 



(57) 



= O (X O1 + X 10 ), 



(58) 
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where 0o is the matrix element < O|0|l > which can be identified with the previous <pi 
in a good approximation. The nondiagonal Hubbard operators X 01 and X 10 are given by 
|0 >< 1| and |1 >< 0|, respectively. The (imaginary) time evolution of the nondiagonal 
operators due to Hi is 

X 01 (r) = e- AET X 01 (0), (59) 

X 10 (r) = e AET X 10 (0). (60) 

Using Eqs.(59) and (60) the displacement Green's function Di(iu n ) can be calculated in a 
straightforward way and yields exactly the previous expression Eq.(53). Taking also Vnl 
into account changes the extremal instanton configurations. However, in the dilute instanton 
approximation, only the action along the extremal path enters. This means that AE acquires 
a different meaning in the effective two-level model but otherwise there is no change. 

In the following it is convenient to use a description based on a Hamiltonian and not on 
a Lagrangian. To this end we consider the following effective Hamiltonian 

H eff = J2 ~ e l a \* a l° + Z! V Mo C <> + C l a la) ~ V Yl 4 C <r 
la la a 

+ ^(X 11 - X 00 ) + gM* 01 + X W ) E *W (61) 

The first three terms in H e ff represent an Anderson Hamiltonian for impurity electrons, 
described by c a} c a , interacting with bath electrons, described by a\ a1 a\ a . As shown in Ref.0 
elimination of the bath electrons allows to represent a general function Go in the effective 
impurity action. The fourth term represents if/, and the last term the interaction between 
impurity electrons and the tunnel modes. 

Neglecting self-energy effects in the tunnel modes the equation of motion for the com- 
posite operator X pq {r)c a {r) is, using H e ff, 

d |x 00 ^ X 01 (-AEc a + ^) . 

^(X^c CT )= dT K 9t) . (62) 

dT { X w (AEc a + d -§?) X 11 ^ 
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with p, q — 0, 1. The equations of motion for the operators c a and ai a are 



<9a 



la 



dr 



= —la\a - ViC a . 



Inserting Eq.(64) into Eq.(63) yields 

rfi .r_. ... d 



J q dr'[5(r - t'K-^j + (J.) - A(r - r')]c a {r') = -gM^ 01 + X 10 ), 



with the function 



A(r-r')=Tj:V l 2 - 

l,n 



Inserting Eq.(65) into Eq.(62) yields 

rP _ .r_. .. . d 



J dr'[5(r - r')( — - + A"(r - r')]x^(r')c a (r') = £ M pq , pW X^ {r)c a {r 
T p'q' 



The only non-vanishing elements of M are 



M 00 ,01 = Mn,oo = Mio.ii = Mn, io = -9<f>o, 



Mio,io = -Mm . m = A£ 



01,01 



We also used the definitions 



A 00 (r - t 1 ) = A n (r - r') = A(r - r'), 



A0 1 ( r - r ')=A( r - r ')e-^(—') ) 



A 10 (r-r') = A(r-r')e 



_ _/n AB(r-r') 



Let us now define the following Matsubara Green's functions 



pq 



G( ™ / ) = -<T(X^(r)c a (r)cUr'))> 

T — T 
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They satisfy the equation of motion 



J dr"G \ r n Tll )G{Jl J + E M P9 , pV G(/_ 9 r/ ) = tf(r - r') < X- >, (74) 
with the definition 

G - 1 ( r m r ,) = S(r - r')(-^ + n) - A«(r - r'). (75) 
Comparing Eqs.(70) and (75) with Eqs.(14) and (15) of Ref0 one finds that 

G ( °° ,) = Go( 11 ,) = G (t-t'). (76) 

T — T T — T 

Furthermore, one can easily verify that Eq.(75) yields 

G ( 10 / ) = G (r-r')e A ^— (77) 
r — t 

G ( 01 / ) = G (r-r')e- A ^'). (78) 

T — T 

Using these results and eliminating the non-diagonal Green's functions one finds that Eq.(74) 
is equivalent with the two equations 

G{ °° ,) = G (r - t') < X 00 > 

T — T 

+g 2 cg f dr"dr'"G Q {r - r")G (r" - r"')e~^"~^G{ w °° ,), (79) 

JO T — T 

G( 11 / ) = G (r-r')<X 11 > 
r — r 

+gYo f dr"dr'"G (r - r")G (r" - r'")e A ^"^G{ w °° ,). (80) 

JO T — T 

Eqs.(79) and (80) can easily be solved in cj-space. One obtains then an explicit expression 
for the wanted Green's function 

G(tuj n )=G( 00 ) + G( n ). (81) 

lUJ n lLO n 

From G one obtains the following expression for the self-energy £/: 
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E/(icj„) 



i - G (iu n )(< x 00 > n( 2 )^ n )+ < > iK 1 )^)) 



(82) 



with 



eT AE? _ j 



(83) 



m 



iu; n - icj m =f AE ' 



where the upper (lower) signs in front of AE refer to the superscripts 1 (2). In the weak- 
coupling limit j < 1, but also in the high-temperature limit AE(3 <C 1, Eqs.(82) and (83) 
reduce exactly to the previous result Eq.(54). The same also holds at zero temperature 
where < X 00 >= 1 and < X 11 >= 0. Eqs.(54) and (82) have, for a finite AE, also the same 
low temperature limit which is a rather astonishing result. 

Taking all together the local Green's function has been obtained in general as a sum of 
three contributions 



Z a is the contribution of the minimum a to the partition function Z. At small couplings 
V has only one minimum at = 0. There exists then only the first term in Eq.(84) with 
Zi/Z = 1 and the perturbation theory can be carried out in powers of as usual. Beyond 
a critical coupling strength V develops two minima. The perturbation theory must then be 
performed around each minimum separately yielding the first two terms in Eq.(58). There 
are in addition tunnelling processes between the two minima causing the third term in 
Eq.(84). The appearance of the prefactors Z a /Z in Eq.(84) is due to the fact that the use 
of the linked cluster theorem in the expansion around a cancels the contribution Z a . The 
definition Eq.(37) of Gi oc implies then the above prefactors. The use of the linked cluster 
theorem in calculating Gj means that there is no prefactor in front of Gj in Eq.(84). In 
general Z = Z 1 + Z 2 + Zj where Zi is the instanton contribution to Gi oc . A convenient way 
to determine Zj is to apply the sum rule to Eq.(84). In practice Zj is small compared to Z 
so that Z\jZ = Z 2 /Z = 1/2 can be used in a very good approximation. 



G L 



oc 



(84) 
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V. NUMERICAL RESULTS 



In this section we discuss the results obtained by solving numerically the set of dynamical 
mean-field equations on the Bethe lattice. < <fi > was determined from Eq.(15) from the 
outset. Then we proceeded in an iterative way, starting with an educated guess for Go- 
The effective potential V((f>) was then calculated from Eq.(21). The calculation of the local 
Green's function Gi oc , given by Eq.(84), amounted to calculating the self-energies S Q and 
£/. To account also for anharmonic effects we determined the stationary solutions of the 
Schrodinger equation for one atom in the potential V: 

^M + 2[E r -V(4>)]M^) = 0. (85) 

r labels the different states with energy E r and eigenfunction tp r (<j>). In the weak-coupling 
case, where V has only one minimum, we obtained the phonon propagator D using the exact 
energies and eigenstates, 

g(i ,„) * Ee - g . jwmi-e-^) m 

^ph r j£- * 'r 



with Q rs — E r — E s and Z p h = ^Z r e~^ Er . Only the first term in Eq.(84) is present and 
the corresponding self-energy is calculated from the diagram in Fig. la) using Eq.(86) for the 
wavy line. As a result anharmonic effects were fully taken into account whereas higher-order 
skeleton diagrams and also Vnl can be neglected as shown in the previous section. 

In the strong-coupling case V exhibits two minima. We treated the two lowest levels as a 
tunnel system and omitted the transitions between them in the phonon propagator Eq.(86). 
Furthermore, the remaining transitions in Eq.(86) can be ascribed either to the left or the 
right well because at least one wave function involved in the transition is localized. We found 
it most convenient to accomplish this splitting of D into left and right parts by dividing the 
matrix element < s\(p\r > into contributions with > and (ft < using also appropriate 
normalization factors. The self-energies T, a ,a = 1,2, were then calculated using again the 
diagram in Fig.l (a) with the split phonon propagator as wavy line. Our procedure becomes 
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exact at large couplings and also provides a very smooth transition to the one-well case. 
Finally Gj was determined from Eq.(82). Since Eq.(82) approaches the simpler expression 
Eq.(52) in several limiting cases we found it sufficient to use Eq.(52) and also to neglect Vnl 
in calculating AE. Having obtained Gi oc we determined a new Go from Eq.(9) and started 
a new cycle until convergency was reached. 

A. Atomic properties 

Fig. 5 shows the self-consistent potential V (in unit of the hopping t) as a function 
of the dimensionless phonon field £ = (fxvo/t 1 / 2 for T = 0.1uj = 0.011 The dotted line 
corresponds to A = 0.5, i.e., the weak-coupling case where V exhibits one minimum. The 
dashed curve refers to A = 1.2 and corresponds to the case where the potential is very flat 
around £ ~ and the double well structure has just been formed. The double well becomes 
more and more pronounced with increasing A as shown by the dot-dashed and solid lines. 
The dotted line in Fig. 6 shows the position of the potential minimum R = ^iLdo/t 1 / 2 as 
a function of A for T = O.Iuq = 0.011 The circles are calculated values which have been 
joined smoothly by the dotted curve. The Figure illustrates how rapidly R increases above 
the threshold value A c ~ 1.15. At large A's R approaches the curve v^A shown in Fig. 6 
as a solid line. The squares in Fig. 7, joined smoothly by a dashed line, represent the 
difference (V max — V m i n )/10 in units of t as a function of A. V max and V m i n are the values 
V(0) and the minimum value of V, respectively. V max — V min approaches at large A the 
curve A/2. The filled circles in Fig. 7 describe the energy difference Q 10 between the two 
lowest eigenvalues, i.e., a phonon frequency. This frequency is at A = equal to the bare 
phonon frequency ujq = 0.11 It softens with increasing A and tends to a very small value 
when the double well has been formed. This "soft" phonon does not produce a structural 
phase transition for A > A c but switches over smoothly to the tunnel mode connecting the 
two potential minima. The corresponding frequency is practical zero on the energy scale of 
Fig. 7. The empty circles in Fig. 7 describe the excitation energy f2 2 o- At A = it is equal 
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to 2uj , softens appreciably with increasing A up to the critical value A c , and then recovers 
approaching u from below at large A's. Fig. 8 depicts three curves for the probability 
function P(£) to find an atom at a distance £ as a function of £. P shows for A = 0.5 
a one-peak structure which broadens and finally splits into two peaks with increasing A. 
Finally, Fig. 9 shows the electro n-phonon correlation function C = —g jt < <pn > versus A 
for T = O.luo = O.Olt. In the weak-coupling limit Co is very small which means that the 
atomic displacements are rather independent from electronic fluctuations and depend only 
on the average electronic density. This changes dramatically for A > A c where the atomic 
displacement depend strongly on the momentary electronic occupation of the atomic site. 

The interpretation of the above results is rather straightforward. In the weak-coupling 
case, for instance at A = 0.5, the atoms are according to their probability distribution 
function located mainly near £ ~ and have only a small mean square displacement from 
this position. From Eq.(9) follows that the charge fluctuations around the average value 
n — 1 are then also small. This is the situation corresponding to the description of a 
metal in the Migdal approximation: The electrons hop very often between sites during one 
atomic oscillation. As a result the atoms oscillate around their unperturbed equilibrium 
positions with only slightly changed frequencies. With increasing A the potential energy of 
the electrons becomes more and more important compared to their kinetic energies and the 
electrons prefer to stay in doubly occupied sites. Since the average occupation of one site 
is forced to be one this means that each site is occupied for long times either with two or 
with zero electrons. The atoms adapt to each of this ionic configuration by relaxing either 
to the left or right potential minimum. If one assumes that the atoms have enough time 
to relax completely to the static equilibrium for each ionic configuration R is given by y/\. 
The solid line in Fig. 6 shows indeed that R approaches y/X at large A's. Using a similar 
consideration V max — V min should approach the curve A/2, which is in agreement with Fig. 
7. Finally, one expects that at large couplings the atoms oscillate more or less with the 
unperturbed frequency around the new equilibrium positions for each ionic configuration. 
This conclusion is also supported by the curve for fl 2 o i n Fig- 7. 
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Our path integral method coincides with the usual Migdal-Eliashberg perturbation ex- 
pansion in the weak-coupling limit. In the strong-coupling case our method resembles in 
several respects treatments based on the Lang-Firsov transformation0'0. One common 
feature is that the atoms acquire new equilibrium positions depending on the momentary 
electronic configuration. In the Lang-Firsov approach the new equilibrium values are given 
by £ = ±y/\. Fig. 6 shows for the Migdal parameter uo^/t = 0.1 that our expansion points 
±R are near to these values for A ^> A c . The perturbation expansion around these points, 
however, is ruled by different expansion parameters: In the Lang-Firsov cas by 1/A and 
in our case by A-uio/t. We also note that the appearance of expansion points different from 
zero signalizes the breakdown of the Migdal-Eliashberg expansion: In case of a double-well 
the information obtained by expansions around the minima cannot retrieved by low-order 
diagrams of an expansion around the local maximum at £ = 0. Our calculations imply that 
the validity of the Migdal-Eliashberg perturbation expansion is not solely governed by the 
condition A ■ u /t <C 1, as often is assumed, but also by A < A c , where A c lies around 1.2 
for a small Migdal parameter Uo/t. Based on the Lang-Firsov transformation the condition 
A < A c for the validity of the Migdal-Eliashberg expansion has also been stressed in Ref.0. 



B. Electronic properties at low temperatures 

In this section we present numerical results for the electronic properties of the Holstein 
model. In the last section we have seen that the crossing from the weak to the strong- 
coupling regime is characterized by a drastic change in the effective atomic potential due to 
the appearance of anharmonic effects and the formation of a double well. These effects also 
cause typical features in the spectral properties of electrons. Fig. 10 shows the evolution 
of the spectral function p(u) = —Im Gi oc (u>)/tt with coupling strength A as a function of 
frequency us for a finite but low temperature T — Q.Iujq — 0.0 It. In case of the weak-coupling 
curve with A = 0.5 p(uS) is somewhat broadened compared to the semicircular density of 
free electronic states. In the metallic state p(0) should always be unrenormalized and equal 



26 



to 1/n at T — 0. Our value for p(0) is slightly smaller than l/V because of the finite 
temperature. For the next lower curve with A = 1.25 the double well has already been 
formed. Nevertheless, p{uj) still consists of a broadened peak centered at u = 0. Atop of 
this broad structure is a narrow component caused also by the tunneling of atoms between 
the potential minima. With increasing A the broad peak for p splits symmetrically into 
two peak centered around ±y/\R with widths comparable to the width of the density of 
bare states. For A < 1.5 the system is still in the metallic state, however, p(0) is much 
smaller than the canonical value 1/tt. The reason for this is that the tunneling contribution 
mainly determines the value of p(0). This contribution, however, is heavily quenched if the 
temperature is larger or comparable with the tunnelling energy AE. The latter becomes 
exponentially small at large A's. For our temperature T = 0.1u the tunnel contribution to 
p{ui) is invisible small in plots like Fig. 10 when the two peaks already separate at around 
A ~ 1.5. The interpretation of the curves in Fig. 10 is simple: For A <~ 1.3 charge 
fluctuations at a site around the average value 1 are small. For values for A between ~ 1.3 
and 1.5 the charge fluctuations become larger and larger and the electronic state at a site 
consists of a coherent superposition of a doubly and an unoccupied state. At A ~ 1.5 the 
coherency is practically lost leading to an incoherent arrangement of doubly and unoccupied 
sites. 

Fig. 11 shows the negative imaginary part of G . For small A's this function resembles 
the density of bare states. With increasing A it develops a strong central component which 
looks like a needle for A > 1.5. At the same time spectral weight is shifted to the upper and 
lower ends of the spectrum which gradually transforms into an upper and lower subband. 
With increasing A these subbands move away from the center and lose spectral weight. Fig. 
12 shows the evolution of — Im S(cj) with A as a function of uo for T = O.Icjo — O.Oli. 
For the values of A on the left panel of this figure — Im X is essentially proportional to the 
density of bare states except at small frequencies where it drops to small values reflecting 
the Fermi liquid ground state. The region of small values decreases with increasing A due to 
the softening of the phonon and the appearance of the tunnel mode. For the larger values 
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for A used in the right panel the spectral weight accumulates in a central component and 
two sidebands which move away from the center. 

The low-frequency behavior of the real part of the self-energy can be characterized by 
an effective mass m* defined by 

dReYi(uj) 



m* jm 



duo 



(87) 



cj=0 

Figure 13 shows m* jm as a function of A for T = O.lcuo = O.Olt. At small A the mass 
enhancement is given by 1 + X/tc. Above A ~ 0.8 m* increases much stronger than linearly 
which announces localization and the transition to the insulating state. The metal-insulator 
transition can also be characterized by the frequency dependence of the ac conductivity cr(fi) 
given by0 

2p 2 tir r r 

a(Q) = — — J deN(e) J dcoA(to, e)A(u + fi, e)(/(w) - f(u + SI)). (88) 

A(u, e) = —Im G(u, e)/7r is the one particle spectral function of the lattice Green's function. 
<y{uo) is shown in Fig. 14 for four different values of A. The large Drude peak seen for A = 1.2 
splits for A = 1.36 into a small peak at uo = and a broader peak at a finite value of to. 
Beyond A ~ 1.4 the Drude peak has virtually vanished and all spectral weight has been 
accumulated in the broad peak caused by transitions between the two subbands. 

Fig. 15 shows the frequency dependence of the spectral function A(u, e) of the lattice 
Green's function at the Fermi surface e = 0. The thick solid line corresponding to A = 1.36 
exhibits the quasiparticle peak at u = and a rather well developed two-peak structure of the 
incoherent background. The thin solid line represents the same spectral function calculated 
within the Migdal-Eliashberg perturbation expansion where the phonon field is expanded 
around = 0. The same parameter has been used as well as the exact eigenstates in the 
potential V. The quasi-particle peak is in this case more pronounced and the incoherent 
background is structureless. Our path integral approach shows that the point = is an 
unstable extremal point in this case and that the two potential minima are stable extremal 
points suitable for expansions. This clearly means that the thin line in Fig. 15 is incorrect 
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and that the Migdal- Eliashberg expansion has broken down in spite of the small Migdal 
parameter 7 = 0.1 used in this figure. For A = 2 both the thick and thin dashed-dotted curves 
exhibit no longer visible quasiparticle peaks at u = 0. The Migdal-Eliashberg expansion, 
however, is unable to reproduce the two-peak structure of the incoherent background but 
instead shows a broad and structureless peak. Though a double well potential has been 
formed already for A = 1.25 the two approaches give similar results in this case, mainly 
because the first excited state is above the barrier and the trajectories of the atom are 
rather delocalized. This means that for a Migdal parameter of 0.1 the Migdal-Eliashberg 
expansion breaks down for A > 1.25. 

C. Electronic properties at high temperatures 

To complete the overview of the electronic properties we show here also results obtained 
in the limit T u . This limit is equivalent to the limit u — > and corresponds to the 
case where the phonons behave like a quasi-static, classical field. Gi oc can then be calculated 
from Eq.(38) and an explicit high-temperature expansion of Eq.(45) and (46) shows that the 
corrections are of the order \u /T. On the other hand, the validity of our general quantum- 
mechanical approach requires \LU /t <C 1 and XT jt <C I. The first inequality follows from 
Eqs.(51) and (52), the second one by expanding Eqs.(45) and (46) at high temperatures 
after taking the limit ujq — > 0. 

The above expansion parameters suggest that the two approaches should yield similar 
results if uoq <C T <C t. To check this we calculated p(uo) for T = 2.0cj = 0.21 The 
results are shown in Figs. 16 and 17. For A = 0.5 the above expansion parameters are 
small and the two corresponding curves indeed are practically coincident. This agrees with 
the observation that in the quantum- mechanical approach Im S(cj = 0) is already so large 
that the quasi-particle peak in the spectral function of the lattice Green's function has 
completely vanished. For A > 1.2 similar double well potentials form in the two cases. In 
the quantum-mechanical treatment the instanton contribution can safely be neglected at 
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these high temperatures and thermally activated processes seem somewhat impeded by the 
level quantization compared to the classical case. As a result the splitting of the band into 
two subbands occurs more slowly in the classical, quasi-static case than in the quantum- 
mechanical case. For A = 4. p{uj) shows in both cases well-splitted bands with a pseudo-gap 
in between where the density is finite but exponentially small. In Fig. 16 the shape of the two 
peaks resembles the bare density of the Bethe lattice whereas in Fig. 17 the shape is more 
Gaussian with well-developed tails into the pseudo-gap. The two Figures illustrate that the 
quantum-mechanical treatment approaches smoothly the classical case with infinite heavy 
atoms at high temperatures and agrees with the latter also quantitatively if all expansion 
parameters are small. 

VI. CONCLUSIONS 

We have studied the Holstein model in infinite dimensions at half-filling in the normal 
state. In contrast to recent semiclassical treatments we consider the atomic displacements 
as a quantized, dynamical field which enables us to recover Fermi liquid properties at low 
temperatures in the metallic state. Our approximation scheme is based on expansions of 
path integrals around extremal paths of the action. This method allows to treat arbitrary 
strengths of the coupling constant A on the same footing. For each stationary path a 
Migdal-Eliashberg-like expansion is generated which converges well both for weak and strong 
couplings. Beyond a critical value of A there is also a time-dependent, extremal instanton 
path associated with the tunneling of atoms between equivalent minima in the effective 
atomic potential. Assuming that the instantons are dilute a perturbation expansion of the 
electronic self-energy in terms of unrenormalized instanton propagators has been carried out 
and shown that for most purposes only the lowest-order, bare term has to be included. 

We found that the usual Migdal-Eliashberg perturbation expansion breaks down around 
A ~ 1.25. This breakdown is caused by the appearance of more than one extremal path of 
the action. Related to this is that the original expansion point at = for the phonon field <f> 
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corresponds no longer to a minimum but to a local maximum of the effective potential. The 
following picture emerges from the calculations: For A < 1.2 the system is in a metallic state 
which can be described by Migdal-Eliashberg theory. Phonon renormalizations are then weak 
and the electronic mass enhancement is linear in A. For A > 1.2 the adiabatic approximation 
breaks down. The system stays mainly in ionic states with either two or zero electrons. Each 
state has its own displacement pattern and instantaneous equilibrium position. The system 
is still metallic due to atomic tunneling processes between the potential minima, however, the 
electronic spectral functions develop a two-peak structure which cannot be obtained in the 
Migdal-Eliashberg theory. Increasing A further the electronic band splits into two subbands, 
the system becomes effectively insulating and behaves like a mixture of unoccupied and 
doubly occupied states. The semiclassical treatment and also the exact solution for vanishing 
hopping show that this "insulating" state is for any finite temperature not truly insulating: 
The density of electronic states between the two subbands is exponentially small but not 
exactly zero. 
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FIG. 1. Skeleton graphs for the electronic self-energy of second (a) and fourth order (b) in g. 
The double solid lines denote the dressed electronic Green's function and the double wavy lines 
denote the full phonon propagator. 
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FIG. 2. Imaginary part of electronic self-energy contributions versus frequency for A 
The solid line refers to the diagram (a) and the dashed line to the diagram (b) of Fig.l. 
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FIG. 3. Imaginary part of electronic self-energy contributions versus frequency for A = 1.375. 
The solid line refers to the diagram (a) and the dashed line to the diagram (b) of Fig.l. 



36 



0.00 



-0.10 



3 

£H -0.20 



-0.30 




A^l.O 
A^l.375 



-0.40 



-2.0 



-1.0 



0.0 

5 



1.0 



2.0 



FIG. 4. Imaginary part of the phonon self-energy T1(ujo) as a function of the dimensionless field 
£ = (Puoo/t 1 / 2 for T = O.lwo = O.Olt. The arrows show the correspondent position of the minima 
of the local effective potential. 
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FIG. 5. Effective potential V of the atom versus the dimensionless phonon field £ = (frujo/t 1 / 2 
for T = O.lwo = 0.01/ and different A's. 
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FIG. 6. Dotted line with circles: Position of the potential minimum R as a function of A, for 
T = O.Icjq = O.Oli. Solid line: Asymptotic A 1 / 2 dependence for vanishing hopping. 
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FIG. 7. Potential barrier V ma x — Vmin (squares) and energy differences between the ground state 
and the first (filled circles) and the second (empty circles) excited states of the effective potential 
as a function of A for T = O.Icjq = O.Olt. 
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FIG. 9. Static electron-phonon correlation function Co = — g /tiffin j versus A for 
T = O.lwo = O.Olt. 
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FIG. 10. Electronic density of states p versus frequency for T = O.Iujq = O.Oli and different A's. 
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FIG. 11. Imaginary part of the "bare" electronic Green's function Go versus frequency for 
T = OAloq = 0.01* and different A's. 
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FIG. 12. Imaginary part of the electronic self-energy X versus frequency for T = O.lwo = O.Olt. 
The chosen values of A extend from the weak- (a) to the strong-coupling (b) regime. 
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FIG. 13. Ratio of the effective and bare electronic mass versus A for T = O.Iujq = O.Olt. 
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FIG. 14. Optical conductivity a in arbitrary units as a function of the frequency for 
T = O.Iuq = 0.01/ and different A's. 
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FIG. 15. One-particle spectral function A at the Fermi surface as a function of the frequency for 
T = O.Iujq = O.Olt. The thin solid and dot-dashed lines are calculated using the Migdal-Eliashberg 
theory, the corresponding thick lines are the present results. 
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FIG. 17. Electronic density of state p versus frequency for T = 2.0ljo = 0.2t and different A's 
in the limit of classical, infinite heavy atoms. 
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